Entropic Tightening of Vibrated Chains 
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We investigate experimentally the distribution of configurations of a ring with an elementary topo- 
logical constraint, a "figure-8" twist. Using vibrated granular chains, which permit controlled prepa- 
ration and direct observation of such a constraint, we show that configurations where one of the 
loops is tight and the second is large are strongly preferred. This agrees with recent predictions 
for equilibrium properties of topologically-constrained polymers. However, the dynamics of the 
tightening process weakly violate detailed balance, a signature of the nonequilibrium nature of this 
system. 
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Topological constraints such as knots and entangle- 
ments constantly form and relax in polymer systems and 
in biomolecules such as DNA . These constraints in- 
crease relaxation time scales and additionally restrict the 
phase space accessible by the macromolecule. Whereas 
the role entanglements play in chain dynamics is well 
appreciated even more basic questions concerning 

effects of topological constraints on static properties such 
as the chain structure remain largely unanswered. Recent 
theoretical studies predict that in equilibrium, a knot- 
ted polymer will generally favor configurations where the 
knot is "tight", i.e., localized to a small region of the 
chain Numerical simulations and scaling analy- 

sis support this prediction pT[ |, but direct experimental 
tests are difficult In this study, we examine this 

interesting prediction experimentally using vertically vi- 
brated granular chains. 

Granular chains consist of spherical beads connected 
by rods, and their backbone enforces the same geomet- 
rical constraints as in a polymer system. A vibrating 
plate supplies the system with energy, balancing the en- 
ergy dissipation due to inelastic collisions experienced by 
beads ||l3[-p^. This polymer system is well suited for 
studying topological constraints as it allows control of 
the chain size and the constraint type, as well as direct 
observations of the chain conformation in contrast with 
traditional polymer systems. Recent studies have suc- 
cessfully used vibrated chains to study diffusive relax- 
ation and spontaneous formation [|l8| of knots. 

We considered the simplest possible topology, a 
"figure-8": a once- twisted ring consisting of two loops, 
separated by a single crossing point, which functions as 
the topological constraint (see Fig. 1). Under appropri- 
ate vibration amplitude, the system is effectively two- 
dimensional and the crossing point hops along the chain 
without ffipping open the figure-8. Surprisingly, despite 
the highly nonequilibrium drive applied to the system, 
we observed strong entropic tightening. The microscopic 
degrees of freedom, the beads, experience periodic drive 
and dissipative collisions with the plate, rods, and other 
beads, as well as frictional forces. Remarkably the macro- 



scopic observable, the loop size, obeys a statistical me- 
chanics. Detailed balance is only weakly violated and 
the empirical loop-size distribution is close to that con- 
jectured on entropic grounds. This provides an opportu- 
nity for studying the role of entropy in nonequilibrium 
systems. We now review the equilibrium result for the 
loop-size distribution. 





FIG. 1. A vibrated ring with a figure-8 twist. 

Tightening of knots in polymers can be understood by 
considering the simplest case, in which electrostatic in- 
teractions are perfectly screened. Suppose a knot in the 
polymer can be considered as located at a point, with sev- 
eral loops of lengths Ni, N2, ... projecting from the knot. 

Ignoring hard-core interactions, each loop has a linear 

1/2 

size of order N^' , so that the number of configurations 
in which the chain returns to itself and forms a loop is 
proportional to e'^^ Yli '^^^j where d = 2 is the spatial 
dimension, and c is a constant. Therefore, assuming all 
microscopic configurations are equally likely, the prob- 
ability of having a given A'^i , N2 , ■ • ■ is maximized when 
Ni K, N and N2, - ■ ■ are all as short as possible. That is, 
the knot spontaneously tightens due to entropic effects. 

A ring with a figure-8 constraint is natural for studying 
this effect as it simply contains two loops of sizes Ni , N2 
with the overall size N — Ni + N2 fixed. The argument 
above predicts a power law divergence in the loop-size 
probability distribution 
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with n = Ni,N2 for 1 <C n -C A^. The exponent a = 
43/16 « 2.7 can be obtamcd even when excluded volume 
interactions are taken into account . Although this 

is obtained from large n, N asymptotics, an enumeration 
for smaller sizes only leads to small corrections to this 
form. 

Our apparatus consists of an anodized aluminum plate, 
driven sinusoidally by an electromechanical vibrator. 
The bead and rod chains consist of hollow nickel-plated 
stainless steel spheres of radius 1.18±0.01 mm connected 
by thin rigid rods of radius 0.26±0.01 mm. The rods con- 
strain both the bending and stretching of the chain. In 
particular, the rods must lie within a cone of a half-angle 
of roughly 25° about the axis of either of the two adjacent 
rods, and the separation b between two adjacent beads 
lies in the range < 5 < 0.94 mm. The chains were con- 
nected end-to-end to form rings, and then twisted with 
a single crossing point thereby forming a figure-8. For 
the experiments reported here, the number of beads in 
the figure-8 ring was between 69 and 219, much larger 
than the tightest possible 8-bead loop. The plate was 
oscillated harmonically at a frequency of 16 Hz. 

The dynamics of the crossing depend on the rms accel- 
eration of the plate, 7 (this dimensionless quantity is in 
units of the gravitational acceleration g). For 7 ^ 1.35, 
the crossing does not move along the chain. For 7 > 1.55, 
the vertical motion of the chain is large enough that the 
number of crossings in the ring is not fixed: a loop of the 
figure-8 can easily flip, untwisting the ring, or creating 
additional crossings. We chose 7 = 1.5, for which flipping 
events remained rare, occurring every roughly 10* oscil- 
lation cycles, while the crossing remained mobile, with 
approximately 50% chance of the loop size changing in a 
1/16-th second cycle. In Fig. 2 we show the distribution 
of change in loop size in one cycle. The acceleration was 
constant and uniform across the plate to better than 1%. 
The 27.2 cm plate diameter, corresponding to approx- 
imately 115 beads, was large enough so that collisions 
with the sloped acrylic wall were rare. 

Digital images of the chain were obtained to determine 
the loop size distribution. Image analysis requires a two 
step procedure involving: (i) monomer recognition, and 
(ii) chain reconstruction. To obtain the monomer po- 
sitions, images of resolution 1000 x 1016 pixels were ac- 
quired. At this resolution, the reflected light from a bead 
appears in the images as a bright spot of about 5x5 pix- 
els, even though the bead has a diameter of just over 8 
pixels. The positions of the beads were determined by 
fitting the intensity pattern generated by each bead to a 
Gaussian, with the peak position taken as the bead posi- 
tion. We estimate the positional accuracy obtained using 
this procedure as 0.05 bead diameters. 
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Change in loop size 

FIG. 2. Distribution of absolute change in loop size over 
1/16-th second, showing the rapidly decreasing probability of 
larger jumps. 

Given these positions, the order of the monomers along 
the chain was determined using an efficient greedy al- 
gorithm which requires only 7V^ operations for an A^- 
bead ring. This algorithm utilizes the aforementioned 
geometrical restrictions on stretching and bending im- 
posed by the rods (given two connected beads, the third 
was searched only in a properly restricted neighborhood) . 
Once the ring was reconstructed, the crossing point was 
identified, thereby determining the loop sizes. 

We first sampled the loop size at a much slower rate of 
0.25 Hz. Starting with two equal-size loops at each run 
for a chain of size N — 149, we obtained the loop-size 
distribution shown as the inset to Fig. 3. This distri- 
bution has a sharp peak located at the smallest possible 
loop size. Hence, the conjectured tightening is indeed ob- 
served. We repeated the experiments using larger beads, 
changing the driving frequency to 13 Hz, and using chains 
of length 49, 69, 99, and 219, as well as using a plate with 
different roughness. Further, many sets of images at ar- 
bitrary phase with respect to the driving were analyzed. 
In all these cases, the same qualitative loop-size distribu- 
tion emerged. There were, however, some quantitative 
differences with the peak height varying by about 20%. 
Since the number of hopping events before flipping is of 
the same order as N'^ (~ 10*), the flipping introduces 
a dependence on the initial conditions, which is more 
pronounced for the longer chains. Even worse, flipping 
events typically occur when the loops are small, leading 
to a reduction of the peak height. Therefore, we used an 
alternate method of determining the loop size distribu- 
tion. 

To find the true peak height, with flipping events re- 
moved, we measured the loop size at a much faster frame 
rate of 16 Hz, and experimentally determined the transi- 
tion probability tij from a loop of size i to a loop of size 
j. To sample tij for all i, the twisted ring was started 
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manually at various equally-spaced loop sizes, and then 
allowed to run for 200 cycles to let the chain equilibrate, 
after which 200 frames were taken to measure Uj. By 
taking 20,000 total frames over 100 separate runs, we 
obtained an accuracy of 10% on individual tij. After 
equilibration, correlations between successive transitions 
were negligible, implying a Markov process. 

Given Markovian dynamics, it is possible to calculate 
the steady state probability, pi, of the loop having size 
i, from the ij.j by performing a Monte Carlo simula- 
tion of the transition process. For a chain of length 
149, we found the distribution shown as the solid line 
in Fig. 3. Compared to that measured directly, this 
histogram is characterized by a sharper peak and con- 
siderable curvature at the center because flipping events 
and dependence on initial conditions were eliminated, re- 
spectively. The dashed line shows the theoretical curve 
Pi oc [i~°'{N — with a = 43/16. The two curves are 
consistent, although the peak is more sharply defined in 
the solid line. Making quantitative statements about the 
relation between the curves would require much larger 
chain lengths N to obtain a sufficient scaling region. Fur- 
ther, the statistical error in pi is most pronounced in the 
center of the histogram where pi is small. Other chain 
lengths have similar histograms, in this case using 0.25 
Hz transition rates. 
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FIG. 3. Loop-size distribution obtained from transition 

rates (solid line), and equilibrium result (dashed line). In- 
set: Loop-size distribution from direct observation, with same 
scale as main figure. 

The transition rates enable us to check detailed bal- 
ance, a sharp test of the nonequilibrium nature of 
the system. For a system in thermal equilibrium, de- 
tailed balance implies a vanishing net flux between 
any two microscopic states, namely Piti_j = Pjtj.i, or 
= In [{piti^i+j)/ (pi+jti+j^i)] = 0. Interestingly, we 
have found that > and /2(i) < 0, namely, short 

jumps tend to tighten the loop more than long jumps. 



This is shown in Fig. 4, where we plot a moving average 
of fi{i) as the solid line, and a moving average of /2(«) 
as the dashed line. The average of is 0.1 ± 0.02 for 
i > 100. 
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FIG. 4. Plot of smoothed /i(i) (solid line), and smoothed 
/2(i) (dashed line). 

As this is an important point, we have made several 
additional checks. We also considered /i(t) -|- /i(z -|- 1) — 
f2{i + 2), which measures the flux around a three point 
neighborhood. This /o-independent quantity was positive, 
demonstrating violation of detailed balance directly from 
the ti_j. Furthermore, we tested detailed balance using 
0.25 Hz transition rates for other chain lengths, and again 
found that short jumps tighten the loop more than long 
jumps. As a final check, we have tested this procedure 
using surrogate data from a simulated process which sat- 
isfies detailed balance to verify that the violation is not 
an artifact of the reconstruction of p from the 

Since it is surprising that an argument based on count- 
ing of states should be relevant in a nonequilibrium sys- 
tem, we now consider a simplified model. We show that, 
even out of equilibrium, there are entropic tightening 
forces, despite violation of detailed balance and possible 
quantitative changes in the size distribution. We consider 
a chain with linear elastic interactions between neighbor- 
ing beads and ignore self-avoidance. We will first consider 
the equilibrium case, with the chain subject to thermal 
forcing and damping, and then generalize to athermal 
drive. Number the beads from 1 to N, and label the po- 
sition of bead i by x{i). Let the crossing occur at beads 
ni, ni -|- 1 and n2,n2 + l. We will compute the forces on 
the crossing, and from this derive an effective dynamics. 

Assume ni < n2, and consider the loop formed by 
beads ni-|-2, m -|-3, n2 — 2. This loop exerts a force on 
bead ni -|- 1 proportional to x{n\ -l- 2) — x{n\ + 1). This 
force tries to tighten the loop. Summing over normal 
modes of the loop, the average force of a loop of length n 
is found to be kT{c—l/n), with c > 0. The fluctuations in 
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this force introduce an effective noise into the dynamics of 
ni. There is also an effective dissipation: if the crossing 
moves to tighten a given loop, the force exerted by that 
loop temporarily increases (both these contributions are 
determined by short distance effects). We then make 
an approximation that ni , n2 can be treated as particles 
of mass M subject to the above forces. The resulting 
dynamics for n = n2 — ni when n << A'^ is, 
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where v is the dissipation parameter and 77 is the noise. 
By the fluctuation-dissipation theorem, noise and dissi- 
pation in (H) are such that n also behaves as a thermal 
particle at temperature T, giving an equilibrium distri- 
bution oc n~'^l'^ . 

Suppose instead that the chain is subject to forcing 
which is athermal with large, non-Gaussian fluctuations 
at short distances, a reasonable assumption given the col- 
lisions with the plate. One can show that the average 
force remains proportional to — 1/n, but the noise in (|^) 
also becomes athermal. The effect of this is most easily 
understood in the overdamped limit of (^). Then, n ex- 
ecutes a biased random walk with a drift of order \jn 
pof , with additional random jumps of varying size due 
to noise. As a result, large jumps are relatively less bi- 
ased compared with small jumps, and detailed balance 
is violated as in the experiment. At large n this gives 
biased diffusion, with the diffusivity determined by the 
mean-square step size. This yields p cx n^", a 7^ djl. 
For smaller n, the probability distribution is determined 
only by the smaller jumps which are more biased, giving 
a sharpening of the peak in the distribution at a width of 
the order of the largest jump size, consistent with obser- 
vations. This behavior is independent of the precise form 
of the noise, as confirmed by our numerical simulations 
of§. 

Finally, we consider the dynamics of constraints in lin- 
ear chains instead of rings. In this case, knots can open 
at the ends of the chain. Consider a linear chain which 
crosses itself at one point, the analogue of the figure-8 
considered here. Let the crossing occur at links ni,n2, 
with < ni < 712 < iV. The points n-i describe a ran- 
dom walk, with boundary conditions n\ < n2, and a bias 
proportional to l/(7i2 — ni). For a > 1, the two walkers 
form a bound state, and the time for the knot to open 
will behave for large N as for a single random walker. 
This contrasts with the behavior found in a larger ac- 
celeration regime, for which experimental measurements 
of knot opening times showed a purely diffusive behavior 
p^ , with negligible entropic interaction between walkers. 
We speculate that the reason for the reduced interaction 
in the larger acceleration regime is that the increased 
drive takes the system further out of equilibrium. 

In conclusion, we have observed a spontaneous tighten- 
ing of topological constraints in vertically vibrated gran- 



ular chains. The significance of this phenomenon is that 
it indicates that the presence of the constraint merely 
reduces the chain length by a fixed amount, rather then 
leading to an extensive size reduction. For equilibrium 
polymers the tightening arises from entropy. Here, due to 
the strongly nonequilibrium drive applied to the system, 
the bead dynamics is athermal and the crossing dynamics 
breaks detailed balance. However, the loop-size distribu- 
tion remains close to equilibrium. This system provides 
further possibilities for experimental examination of the 
role of entropic forces in nonequilibrium statistical me- 
chanics. It is possible to probe fluctuation-dissipation re- 
lations by a quantitative comparison between the forces 
on the crossing point and the velocity fluctuations in the 
beads, namely the granular temperature. 
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